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Abstract 

We investigate the possibility that a heavy scalar field, whose mass exceeds the Hubble scale 
during inflation, could leave non-negligible signatures in the Cosmic Microwave Background (CMB) 
temperature anisotropy power spectrum through the parametric resonance between its background 
oscillations and the inflaton fluctuations. By assuming the heavy scalar field couples with the 
inflaton derivatively, we show that the resonance can be efficient without spoiling the slow-roll 
inflation. The primordial power spectrum modulated by the resonance has a sharp peak at a 
specific scale and could be an origin of the anomalies observed in the angular power spectrum of 
the CMB. In some values of parameters, the modulated spectrum can fit the observed data better 
than the simple power-law power spectrum, though the resultant improvement of the fit is not 
large enough and hence other observations such as non-Gaussianity are necessary to confirm that 
the CMB anomalies are originated from the resonant effect of the heavy scalar field. The resonant 
signatures can provide an opportunity to observe heavy degrees of freedom during inflation and 
improve our understanding of physics behind inflation. 



I. INTRODUCTION 



Cosmic inflation is the standard paradigm that provides the initial conditions for struc- 
ture formation and the anisotropics in cosmic microwave background (CMB) as well as the 
global properties of the spacetime [IH3]. In the inflationary scenario, the accelerated expan- 
sion stretches quantum fluctuations on microscopic scales to cosmological scales, providing 
the seed for macroscopic observables as the anisotropies in the CMB. The wavelengths of 
the fluctuations are extremely short in the earliest epoch. Thus the cosmological observa- 
tions provide a window into short-distance physics which is beyond the reach of terrestrial 
experiments. 

The cosmological fluctuations, which have been generated quantum mechanically, are 
statistical in nature. In the simplest single-field slow-roll inflation models, they are approx- 
imately Gaussian-distributed and their power spectrum is nearly scale-invariant. Then we 
usually characterize them over the observable range of scales in terms of a power-law-type 
power spectrum, which is parametrized simply by two or three parameters: the amplitude, 
the spectral index, and the running. Though the two- (or three-) parameter spectrum con- 
sistently explains the observed CMB anisotropies [4j, it could miss valuable information on 
physics behind inflation. Recent high-resolution CMB data already implies the presence of 
fine features in the primordial power spectrum. Indeed, several groups, including one of the 
present authors (JY), reported statistically significant discrepancy between the prediction 
from a power-law primordial power spectrum and the CMB data [SJE]- Using non-parametric 
reconstruction methods, they have found large anomalies in the reconstructed power spec- 
trum, which are localized around wavenumbers k ~ 0.003 Mpc -1 and k ~ 0.009 Mpc~ x . On 
the other hand, a number of effects that cause deviations from a power-law power spectrum 
have been investigated in literatures, including trans-Planckian effects [T— 9J , a burst of par- 
ticle production [TOl [TT] , temporal violation of slow-roll approximation [T2HPT] , turns in the 
inflationary trajectory [T8H20] . a sharp water field transition [21], and a sudden change of 
sound velocity [25] . They modify evolution of the inflaton fluctuations in their own way and 
leave their characteristic signatures on the primordial power spectrum. Thus fine features 
in the power spectrum could contain rich information such as detailed structures of inflaton 
Lagrangian or existence of other degrees of freedom. 

In this paper, we consider an effect of coherent oscillations of a heavy scalar field whose 
mass exceeds the Hubble scale, m ^> H. Heavy scalar fields are ubiquitous in models of 
inflation embedded in supergravity and string theory. They appear as moduli fields, Kaluza- 
Klein modes, the scalar supersymmetric partner of inflaton, or others. In usual treatment, 
the dynamics of a heavy scalar field is neglected assuming that it is stuck to their potential 
minima during inflation because its excitation decays quickly [26J. The instantaneously- 
excited oscillations, however, can be important when we measure the primordial power 
spectrum with high resolution. An impact of the excitation has already been discussed 
in Ref. [IS] for a hybrid-inflation model. More recently, Refs. J2DH22] have discussed that 
oscillations excited by a sharp turn in multi-dimensional potential leave a ringing signature in 
the primordial power spectrum. In both cases, the signatures arise because the evolution of 
the fluctuations becomes non-adiabatic [7] by a sudden energy transfer between the inflaton 
and a heavy scalar field through their couplings in the potential. In addition, it has been 
discussed in Ref. [23] that a ringing feature in the power spectrum is induced through the 
gravitational couplings without considering direct couplings between inflaton and a heavy 
scalar field. 
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In this paper, instead, we point out that a resonant enhancement of the fluctuations 
efficiently occurs deep in the horizon, k/a ~ m 3> H, through derivative couplings with a 
heavy scalar field. The derivative couplings are allowed even if a shift symmetry is imposed 
to ensure the flatness of the inflaton potential. Then there is no reason why they are absent 
in the action from the effective-field-theory point of view [27J ES]- Though the derivative 
couplings are usually irrelevant at low energy scale, they can play an important role on the 
evolution of the inflaton fluctuations in the resonance epoch. This is because the inflaton 
fluctuations and the heavy scalar field rapidly oscillate in the resonance epoch. Hence, effects 
of the derivative couplings on the evolution of the inflaton fluctuations are large there while 
their effects on the background evolution is relatively small because the derivative of the 
background inflaton field is slow-roll suppressed. In the following sections, we estimate the 
enhancement of the primordial power spectrum by the resonance assuming non-derivative 
couplings are sufficiently suppressed. In contrast to other effects as the slow-roll violation 
[29] . the feature induced by the resonance can be sharp and large even in the case that 
adiabaticity is mildly violated because the resonance coherently accumulates small effects. 
The instantaneously excited oscillations do not affect the slow-roll background evolution 
much in this case unless the heavy scalar field dominates the energy density because the 
flatness of the inflaton potential is ensured even during the oscillations. 

The organization of this paper is as follows. In ^TTj we present our model to realize 
an efficient enhancement of the fluctuations and discuss conditions on the model required 



for the slow-roll inflation. In £111, we estimate the enhancement of the fluctuations by the 
parametric resonance and discuss consistency of our model with the anomalies observed in 



the CMB spectrum. Finally, we provide our summary of this paper in £IV 



II. THE MODEL 

In this section, we introduce our model and discuss conditions on model parameters 
required for the slow-roll inflation. 

A. An inflationary model with a heavy scalar field 

We consider a model with a heavy scalar field with mass m 3> H which couples to the 
inflaton field in the following manner, 



m 



-(d<j>Y + V{<j>) + -(dxY + —X 2 + K n + K d 



(1) 



with 



(2) 



and 



(3) 
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where is the inflaton field and x is the heavy scalar field. The potential V(<f>) is assumed 
to be sufficiently flat; 



e v < 1, It/vI < 1, (4) 

where 

Ml fV'V V" 
C ^lt{v) • Vv = Ml T , (5) 

are the slow-roll parameters. Here, M p = 2.4 x 10 18 GeV is the reduced Planck mass. 
The heavy scalar field x is assumed to be subdominant; 

fx « !» (6) 

where 

, = Px „ X 2 + m 2 X 2 

Jx p QM 2 H 2 ' 1 ; 

is the fraction of its energy density to the total one. We assume that x decays with a rate 
T, which satisfies H C T <C m. 

The derivative couplings K n and K d are expected to appear in the action if we consider 
the action ([I]) as the leading terms in a generic effective field theory EE]- Though 
other terms are also allowed in general, the derivative couplings K n and Kd provide the 
most general couplings between the inflaton and the heavy scalar field at the leading order 
in 1/A n and 1/A^ in a model with a parity symmetry, — > —(f), and a shift symmetry, 
— > + c. Moreover, the couplings Kd cannot be forbidden while the couplings K n can be 
suppressed by imposing a shift symmetry on x- Higher-order terms in 1/A n and 1/A^ are 
also expected to appear in general, though we have suppressed them in Eq. Q. To ensure 
that contributions from these terms can be safely neglected, we assume hereafter that the 
background fields satisfy the following conditions, 

X < A n , 0, x < AJ. (8) 

We have also suppressed terms (<90) 4 and (dx) 4 , because they have little effect on our analysis 
under the conditions above. In general, non-derivative couplings as could appear since 

the potential breaks the shift symmetry. However, we assume the non- derivative couplings 
are sufficiently suppressed in the following analysis for simplicity. Even if we include non- 
derivative couplings, they modify only the background evolution and do not directly affect 
the evolution of fluctuations much. 

We can also find the couplings K n and Kd in specific models of inflation. For example, 
the couplings K n naturally appear if the inflaton is a pseudo-Nambu Goldstone boson, 
where the heavy scalar field is provided by the symmetry breaking field. Moreover, they 
can be found in the Einstein-frame action if a model contains additional scalar degrees of 
freedom with a non-minimal coupling to curvature, and also in the supergravity action with 
higher-order terms in the Kahler potential. In the latter case, the scalar supersymmetric 
partner of inflaton can be a candidate for the heavy scalar field, which usually has large 
mass corrections from the Planck suppressed terms in the F-term potential since it is not 
protected by a shift symmetry. On the other hand, the couplings Kd can be found in models 
such as brane inflation. In brane inflation, for example, we can see that they appear in the 
DBI action by expanding its square root with the heavy scalar field being provided by KK 
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modes or brane coordinates other than the inflaton. In the case of brane inflation models, 
where all the higher-order terms neglected in the action QlJ) are known, we may not need to 
impose the conditions ^ though instead the evolution equation should be solved including 
their effects. In addition, they have non-derivative couplings and hence we should make the 
analysis including them. These cases will be discussed in another paper. 



B. The background evolution 

The evolution equations for the background fields can be written as 



tt x + (3H + T)n x + m 2 x + 



A 



n / 2 



2A T 



0. 



where 7r^ and n x are conjugate momenta of the scalar fields, 



7I> 



1 + X n ^- n + (A* + A d2 ) ^ 



7T V 



1 + (Adi + A 



J2) 



2A^ 



X 



■, l 



= K 2I cj> . 

Here we have introduced a notation that 



K = \{d4>f + \{d X f + K n + K d , 



(9) 
(10) 

(11) 
(12) 

(13) 
(14) 

(15) 



and 4>^ = 0, 4>^ = x f° r brevity. 1 Kjj represents K differentiated by X IJ = —d4> I -d4> J /2. 
A dissipation term, Tn x , has been added in Eq. (10) to incorporate the effect of the decay 
of X- 



We can obtain the solutions for the background evolution as usual, 

Xif) - Xoe" r *cos(mt), 



(slow-roll solution) 



(16) 
(17) 



if the conditions (|8j) are satisfied at the onset of the oscillations, t = 0. We discuss consistency 
of these solutions in some details here. 



Before proceeding to a discussion on the consistency of the solutions (16) and (17), we 



discuss whether inflation occurs or not during the oscillations. Inflation is realized if the 
slow-variation parameter is small; 



H 
IP 



<C 1. 



(18) 



1 We use the summation convention for repeated indices throughout the paper and parenthesis for compo- 
nents in the field space. 
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Using the Friedman equation, H = — KjjX ij /M^, the equation (18) can be rewritten as 

KjjX ij 



en 



M£H 2 
1 



1 



2 \ M p H 



X 



2 \M p H 



(19) 



where we have neglected the higher-order terms in 1/A„ and 1/A^. Again, using the Fried- 



man equation, ?>M 2 H 

as 



p ~ V, and the solution (17), we can further rewrite Eq. (19) 



e H ^ ey + -/ x sin 2 (mt), 



(20) 



where we have neglected the higher-order terms in ey and f x . Hence, the conditions 
(J6|, and ^ are sufficient for realizing inflation. Note that another slow- variation parameter 
rjH = en I Hen is not so small during the oscillations due to the second term. 

In order that the slow-roll solution (16) be consistent, the first term should be smaller 
than the others in Eq. -k^/SHir^ <C I. Substituting the solution (16), we obtain 



n<t> 



V" <p H 
V 7 !! ~ IP 

-Vv + ch, 



(21) 



which shows that the solution ( 16 ) can be used consistently under the conditions Q, (|6]), and 
It is noteworthy that we cannot replace tt^ by in deriving Eq. (21 ). The differentiation 



of x i n Eq. (11) induces a large factor m, hence 0/3if0 can be large, 



3H<f) 



O 



H A„ 



O 



m x 2 



(22) 



Finally, tt x can be replaced by x i n solving Eq. (10) because both and x induce the factor 



m. Hence the couplings do not affect the solution (17). On the other hand, the last term 



in Eq. (10) induces an approximately constant term in x as 



xit) 



2m 2 A r 



+ Xo + 



2m 2 A r 



-rt 



cos(mt). 



(23) 



Though the constant term does not spoil the resonance, it should satisfy the condition (|8 
This condition can be satisfied if m is sufficiently large, 



m 
H 



K 

M r , 



(24) 



where eh^ is the contribution of the inflaton to en- The constant term also contributes to 
the potential energy. However, the induced potential energy is always sub dominant if the 
condition (24) is satisfied. In the case that the constant term is comparable to the initial 



amplitude Xoj f ne amplitude of the oscillations is provided as Eq.(23) instead of Xo 
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In summary, the oscillations of the heavy scalar field x do not spoil the slow-roll inflation 
provided that the conditions Q, ([6]), and ^ are satisfied. In the next section, we investigate 
the evolution of the inflaton fluctuations in this background and show that a parametric 
amplification of these fluctuations takes place through the resonance with the oscillations of 
X- 



III. PARAMETRIC RESONANCE WITH THE OSCILLATIONS OF THE HEAVY 
SCALAR FIELD 

In the previous section, we saw that the oscillations of the heavy scalar field do not much 
affect the background evolution provided that the conditions Q, ([6j, and ^ are satisfied. 
In this section, we show that the fluctuations in the inflaton field can be enhanced through 
the parametric resonance with the oscillations of \ even in this case. 



A. Evolution equation for the inflaton fluctuations 

First, we derive an evolution equation for the inflaton fluctuations. The heavy scalar field 
X oscillates with a frequency m ^> H, then the resonance occurs on scales much smaller than 
the horizon scale. Hence, we can neglect contributions from the metric fluctuations during 
the resonance because they are much smaller than those in the scalar fields on subhorizon 
scales as explicitly shown in Appendix |Aj Furthermore, the fluctuations in the heavy scalar 
field can be neglected because they do not contribute to the resonance. To eliminate the 
gauge degrees of freedom in the metric fluctuations, we employ the flat gauge, where the 
spatial metric becomes a 2 5ij. 

Neglecting the contributions from the fluctuations in the metric and the heavy scalar 
field, the second-order action in the inflaton fluctuations, ip, can be written as, 

S 2 ~ J dtd 3 x y [{2K lKM X KL + K n )cp 2 - K n {Vy) 2 /a 2 - K ltl cp 2 ] (25) 

* f dtd 3 * | [0 2 - c 2 s (V^/a 2 ] , (26) 

where 



a\2K lKM X 

„3 



KL 



c e = 



2K 1KM X KL + K n 
1 + (Adi 



2A d2 ) - () 



2A^ 



(27) 
(28) 

(29) 

(30) 



Here, we have neglected the potential term, K\^p 2 ~ —3f]vH 2 ip 2 , which is much smaller 
than the term (V(f ) 2 /a 2 during the resonance. Note that non- derivative couplings, if any, 
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can be similarly neglected unless the slow-roll conditions Q are violated. The action (26) 
leads to the following evolution equation in Fourier space for the inflaton fluctuations, 



Vk + 



Vk = 0, 



(31) 



where v = z^tp. Neglecting the higher-order terms in x/A n and the quantities c 2 and 

'z^/z^ are estimated to be 



c^l + (A dl -2A rf2 )^sin 2 (mt), 



(32) 



Zs 



-(3 - 2t H )H 2 + 3mH 



2 «2 

Xn^- sin (mi) + (X dl + 2X d2 ) T ^-^ sin(2mt) 



+ m" 



m 



2 ~ 2 

A n ^- cos(mt) + {X d i + 2A d2 )^^ cos(2mt) 

H 



2 "2 

An oT" cos ( mt ) + (-^i + 2A d2 )^r|^ cos(2mt) + O 



(33) 



where Xo = Xo e r *- Substituting Eq. (32) and Eq. (33) into Eq. (31), the evolution equation 
can be rewritten in the form, 



Vk + vnf 



j _ S^. cos(mt) — 2q d cos(2mt) 

am J 2 



v k = 0, 



(34) 



where 



A, 



Xo 
A„ 



"(A d i - 2A 



r/2j 



2 " 2 

m Xo 



fc 

am 



+ (X d i + 2A 



</2y 



2 ~ 2 

m Xo 
4A^ ' 



(35) 
(36) 



Hereafter, we consider the limit where either q n or q d is much larger than the other for 



simplicity. In this case, the evolution equation (34) can be written in the form of the 
Mathieu equation, 

d 2 v 

+ [A k - 2q cos(2z)] v k = 0, (37) 

where (z,Ak,q) = (mt/2,(2k/am) 2 ,q n ) for the i^-dominated case and (z,Ak,q) = 
(mt , (k j am) 2 , q d ) for the -fQ-dominated case. Here, the value of q is smaller than unity 
because x/A n "C 1 and mx/A^ ~ x/A^ 1 as assumed in the previous section, and then 
the adiabaticity is mildly violated. Hence the resonance is narrow, which occurs in a nar- 
row instability band, \Aj. — 1| < q. The parametric resonance occurs when the modes are 
redshifted to this instability band; 



_ q^\ m k / qn\ m 
2 J 2 a v 2 ) 2 ' 



(38) 
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for the i^-dominated case and 



Qd 
2 



k 

m < — < 
a 



1 + 



Qd 



m, 



(39) 



for the -fQ-dominated case where 



Qd = (A 



di 



6A 



d2) 



1 "1 

™ Xo 
8A 4 , ' 



(40) 



In both cases, the resonance occurs in a deep subhorizon regime since the mass scale m 
is much larger than the Hubble scale, m ^> H. Though we consider only the case that 
one dominates the other, the resonance is expected to be efficient even when both of the 
couplings are non-negligible since the values of their oscillation frequencies are rationally 
related [30J. 

After the oscillations damp out, we can use the standard evolution equation for the 
inflaton field in a single-field inflation model, which can be written as 



Vk + 



MH 2 



v k = 0, 



(41) 



where 



M = -(3-2e H ) + 
3, 

= -(3-2e H ) + 



a 3 M^H 2 dt 



a 3 7r! 



H 



1 2 , VH 

4 Vh + 2H 

Yl 



(42) 
(43) 



in our notation. Here, we have added the contributions from the metric fluctuations and 
the potential term, which induce the global tilt of the spectrum. 



B. Parametric amplification of the curvature perturbations 

We estimate the amplification of the power spectrum for the comoving curvature per- 
turbations relative to the standard one. The comoving curvature perturbations, (, are 
proportional to the inflaton fluctuations in the flat gauge, (p at linear order as £ = Hip /(f) 
after the oscillations damp out. Hence, the ratio of the modulated power spectrum to the 
unmodulated one, A, can be written by the mode function for the inflaton fluctuations as, 



A = lim 


Vk 


2 


= lim 


Vk 


t— >oo 


Vk 




t— >oo 


Vk 



(44) 



where we have denoted the modulated quantity by a symbol with tilde. We present the 
results for the fQ-dominated case here. The analysis for the i^ n -dominated case can be 
made in a similar way. 

The spectral shape of the modulated power spectrum can be roughly understood from 
Fig. [TJ The spectrum has a peak around the modes which had the mass scale ~ m at the 
onset of the oscillations and the peak width for smaller and larger wavenumbers, and Ax, 
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FIG. 1: Time evolution of the infiaton fluctuations. The fluctuations on observable scales can 
be enhanced through the parametric resonance when they are deep in the horizon. The modes 
with comoving scales represented by the thick line are enhanced by the resonance, while those 
represented by the thin lines are not because they never cross the resonance band during the 
oscillations of the heavy scalar field. Hence, the resonance induces a peak in the power spectrum. 
Here, the variables with subscript indicate those evaluated at the onset of the oscillations. 



are respectively determined by the width of the resonance band go an d the duration of the 
oscillations H/Y. From Fig. [TJ we can see that the modes whose wavenumbers are around 
k p / do = m(l + go/2) are most amplified while the modes whose wavenumbers are less than 
ks/cLQ = m(l — go/2) are hardly amplified since they were already outside of the resonance 
band at the onset of the oscillations. Hence, the peak width for smaller wavenumbers is 
roughly given by, 



Ac 



go- 



(45) 



On the other hand, the modes whose wavenumbers are larger than ki/a = e H ' r m are also 
hardly amplified because they cross the resonance band after the oscillations damped out. 
Hence, the peak width for larger wavenumbers is roughly given by, 



A; 



H 

Y' 



(46) 



Moreover, the peak amplitude can be roughly estimated by using the Floquet exponent 
\i = g/2, as, 



Vk oc exp l / fidz 
where the exponent is approximately given by, 



(47) 



mq 

4r 2 

mgg 



AH 

W mm.g , r 



for g > H/Y, 

for g < H/Y, 
H s 



(48) 
(49) 
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Here, the integration range is determined by the condition \Ak — 1| < q- 

In Fig. [2j we have shown the amplification of the power spectrum, A, estimated numer- 
ically assuming that the inflaton fluctuations were in the Bunch Davis vacuum before the 
oscillations start. Each panel shows the spectra for different values of a parameter while the 
others are fixed. The resultant spectra exhibit the expected dependence on the parameters 
and rapidly oscillate with frequencies of the order of m/H. In the i^-dominated case, the 
spectrum behaves similarly. 




FIG. 2: The amplification of the power spectrum for different values of the parameters, m (top left 
panel), V (top right panel), and qo (bottom panel) for the ^-dominated case. The spectra behave 
similarly in the i^ n -dominated case. 

Since the width and the amplitude depend on the model parameters in different ways, we 
can extract information on the heavy scalar field from the features in the primordial power 
spectrum. For example, we can read values of the model parameters in the ^-dominated 
case from the data as, 

r ~ 10,3 &v Got4v) (£tr • (5 ° a) 
^ 10,5 1/2 © 1/4 (^r< <«*> 



ii 



assuming A^i and \d2 being 0(1). Here, we have made a rough estimation of the width as 
As ~ % and Al ~ H/Y and the amplification at the peak scale A p has been estimated as 
Eq. (49). We have also introduced A min = min(As, A L ) for brevity. 



C. Comparison with observational data of CMB 

Having clarified the dependence of the modulated power spectrum on our model parame- 
ters, we can in principle fix or constrain their values by confronting them with observational 
data. In practice, however, it is difficult to find best-fit parameters for this kind of highly 
oscillatory spectrum because there are many minima in the likelihood surface and then the 
Markov Chain Monte Carlo chains do not converge [321 E3] . Moreover, the parameter search 
is time-consuming because we do not have the exact analytic expression of the modulated 
power spectrum. Therefore, instead of searching the best-fit parameters, we just show one 
possible values of the parameters that provide better fit of the observational CMB data than 
the power-law primordial power spectrum. 

In Fig. |3j we have shown the angular power spectrum of the CMB for the modulated 
power spectrum with a feature at k ~ 0.003 Mpc -1 , which corresponds to multipole I ~ 40. 
As possible values of the parameters, we have used m/H = 10 4 , Y/H = 5 x 10 2 , and 
go = 0.1, while the other cosmological parameters have been kept fixed to the WMAP 7 bet- 
fit values [J|. Using these values of the parameters, we have found that the x 2 -value impoves 
2.5 in comparison to the simple power-law spectrum. Hence, the modulated spectrum can 
provide a better fit of the data, though the improvement is not enough to confirm that 
the CMB anomalies are originated from the resonant effect of the heavy scalar field since 
we have introduced three additional parameters. On the other hand, for the feature at 
k ~ 0.009 Mpc -1 , it is difficult to obtain significant improvement. The angular power 
spectrum of the CMB is provided by a convolution of the primordial power spectrum with a 
radiative transfer function, which is written in terms of the spherical Bessel function. Since 
the radiative transfer function for multipole I has an oscillatory tail for larger wavenumbers, 
k ■ (1.4 x 10 4 Mpc) > /, the peak induced in the angular power spectrum has a relatively 
large width even if the primordial power spectrum has a sharp peak. Though the observed 
sharp and high peak could be realized if the modulated power spectrum has an oscillatory 
component coherent with that in the radiative transfer function, this is not the case for the 
spectra obtained here. 

Though the evidence found here is not statistically significant, this will give us some 
insight into the possibility to find some signatures of scalar fields much heavier than the 
Hubble scale. If we could found such signatures, they will provide a hint for physics behind 
inflation. 

IV. SUMMARY AND DISCUSSION 

In this paper, we have discussed the possibility that a heavy scalar field, whose mass ex- 
ceeds the Hubble scale, m 3> H, could leave non-negligible signatures in the CMB spectrum 
through parametric resonance between its background oscillations and the inflaton fluctua- 
tions. The resonance could be efficient without spoiling the slow-roll inflation if the heavy 
scalar field couples with the inflaton derivatively; the feature induced by the resonance can 
be sharp and large even in the case that adiabaticity is only mildly violated because the 
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FIG. 3: The angular power spectrum of the CMB for the modulated power spectrum. 

resonance coherently accumulates small effects. In the analysis, we have assumed that the os- 
cillations of the heavy scalar field are instantaneously excited at e- folds N p ~ ln(m/H) + iV*, 
where N* is the number of e-folds from the end of inflation when the peak modes crossed 
the horizon. This assumption will be appropriate if inflation begins at the e-folds N p and 
the heavy scalar field is initially displaced from its minimum, which is natural in the cases 
that the minima differ before and during inflation, or inflation occurs after a tunneling from 
a neighboring minimum [34H38] . for example. As another possibility, the oscillations could 
be dynamically excited when the heavy scalar field becomes momentarily light/tachyonic 
or the slow-roll condition is temporary violated. In the latter case, effects of the slow-roll 
violation on the resonance should be taken into account to consider large excitations of the 
heavy scalar field [39]. We will address this issue using a specific model in another paper. 

We have also estimated the goodness-of-fit of our model with the anomalies observed in 
the CMB spectrum. The resultant improvement of the fit has not been large enough to 
confirm that the CMB anomalies are originated from the resonance, though a systematic 
analysis has not been performed for technical difficulties. To test our model further from 
observations, non-Gaussianity in the CMB anisotropics will be helpful. Since oscillatory 
components are induced in the interactions, the higher-point correlation functions are also 
enhanced at specific scales by the resonance [4*0Tl4*3"] . We will estimate the amplitude and 
the shape of this non-Gaussianity in an upcoming paper. 
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Appendix A 



Here, we examine the two assumptions made in § |III[ i) the metric fluctuations can be 
neglected during the parametric resonance, where k/a ~ m ^> H and ii) the fluctuations in 
the heavy scalar field, Xi can be neglected. 

To compute the second-order action, we work in the ADM formalism, 



ds 2 = -N dt + 7^ (da; 1 + N l dt)(dx J + N J dt). 



(Al) 



Using the metric (Al), the action for the gravity and the scalar fields can be written as 

'Ml 



S = / dtdx'N^ 



M 2 

_p i?(3) + -^(KyK* - K 2 ) + P(X IJ , d> T ) 



(A2) 



where is the three dimensional Ricci curvature and is the extrinsic curvature of the 
time slice, 

Kij^T^dij-Niti-N^), (A3) 

Here, | denotes the covariant derivative for the spatial metric 7^. We have written the terms 
for the scal ar fie lds simply as P(X IJ , cp 1 ) because a detailed form is not important here. For 
the metric (Al), X IJ is written as 



X 



1 j 



(A4) 



The Friedman equation and the evolution equations for the background fields can be written 
as, 



3M 2 H 2 = p, 

7Tj + 3H7V! -P Z = 0. 



(A5) 
(A6) 



where we have introduced the energy density of the scalar fields, p = 2PjjX J — P, and the 
conjugate momenta of the scalar fields, tvj = Pij(f> . 

We estimate the contributions from the metric fluctuations and the fluctuations in \ by 



expanding the action (A2) explicitly up to second order in the fluctuations. To eliminate 



the gauge degrees of freedom in the metric fluctuations, we adopt the flat gauge condition, 

7y = a 2 5ij, (A7) 
where we have neglected the vector and tensor perturbations. Expanding N and N l as 



N = 1 + a, 

Ni = ft/3, 



(A8) 
(A9) 



second-order terms dependent on these metric fluctuations in the action (A2) can be written 

as 



dtda; 3 a 3 



( PlJ X IJ - p)a 2 - aputfcp 3 - ap^ 1 + ( vrV - 2M 2 Ha) 



d 2 (5 



, (A10) 
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where (p 1 is the fluctuations in (f) 1 . This action leads to the following constraints at the first 
order of the fluctuations, 



a 



7Tj i 



2M 2 H 



d 2 /3 PIJ X 



ij 



P 



M 2 H 



a + 



Pu 
2M 2 H 



<p'<p J - 



Pi 



2M 2 H 



(All) 
(A12) 



fluctuations can be written as 



Substituting Eq. (|A11|) and Eq. (|A12|) into Eq. (|A10|) the contributions from metric 

(A13) 



2 L 



where 



M 



(g) 
ij 



2 V 



( PklX kl 



P\ 7Tj 7Tj 



P 



3 / M p pp\ 7Tj) 



M p H M p H 2 



d (a 



p yMN TTj 

PMN,KLX MpH MpH ~ M P P KL,(I\ 



K 



P 

M p H 



MpH 2a 3 M p 2 H 2 dt \H 



(A14) 



/ X KL 

\Mm< 



+ 



d fa 



a 3 M 2 H 2 dt \H 



(W J ) 



«V(g) _ 9 p x-LAf ^ 

jv n = AriK.LM-^ - - „ — — 



MpH MpH' 



1 1, j = IK,LM 

which coincide with terms obtained in Refs. 



7T/7TJ 

(A15) 
(A16) 



5\. Here, Pij,kl and Pij^k represent 
another derivative of Pu with respect to X KL and <p K respectively. These terms con- 
tribute to the equation of the motion for the fluctuations as Af^H(p J and [Xi'fjH 2 — 



a 3 (a 3 
to be 



Af'f] H)' /2}ip J . At the leading order in e, f x , (2) /A„, and 4> /A 2 ,, Aif) is estimated 



M$ ~ 0(e), M 



(g) 



M„ 



0{e* 12 ) + ?0(e 1/2 /y 2 ), M 



12 An" v ^ ' ' H 



(g) 

22 



m 
H 



0{f x 



(A17) 



On the other hand, jifjf is estimated to be 



Mf ~ AT 2 ? ~ |jO( e ), ~ g [0(e) + 0(f x )) , 



A 4 



(A18) 



and (A/jjp)' ~ mj\f u \ Hence, the contributions from the metric fluctuations are much 
smaller than the term d 2 ip ! /a 2 , which is estimated to be d 2 ^ 1 /a 2 ~ rn^ip 1 at k/a ~ m. 
On the other hand, second-order terms dependent only on the fluctuations in the scalar 
fields are given by 



dtd 3 x 



{2P IK)LJ X KL + P/j)^ V " Pijd i cp I d i ip J /a 2 - Pi,j<p T <p J + 2P K , u '4> K ip 1 Cp J 

(A19) 
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The fluctuations in the heavy scalar field, <fi( 2 \ oscillate with a frequency larger than m 
because of the mass term. Therefore, cp^ does not lead to the parametric resonance in 
contrast to cp^\ Moreover, the mass term also suppresses the amplitude of cp^ at k/a ~ m 
compared to that of <p^ l \ Hence, we can safely discard the fluctuations in the heavy scalar 
field, ^ (2) - 

Neglecting the contributions from the fluctuations in the metric and the heavy scalar field, 
the second-order action reduces to a rather simple form, 



S 2 = J dtd 3 x j [{2P lKM X KL + P n )y? 2 - P n (V^) 2 /a 2 - P^ 2 ] . 



(A20) 



The potential term can be rewritten as P\^ 2 — —3r]yH 2 ip 2 . Hence we can also neglect this 
term. Using the form (IT]) for P, we obtain the action (26) in {111 
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